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Abstract. We investigate crossing path probabilities for two agents that move randomly in 
a bounded region of the plane or on a sphere (denoted R). At each discrete time-step the agents 
move, independently, fixed distances di and d^ at angles that are uniformly distributed in (0,2-7r). 
If Ft is large enough and the initial positions of the agents are uniformly distributed in Ft, then 
the probability of paths crossing at the first time-step is close to 2d\dil \isA\R\), where A[R] is 
the area of R. Simulations suggest that the long-run rate at which paths cross is also close to 
1d\di I '(tt A[i?]) (despite marked departures from uniformity and independence conditions needed for 
such a conclusion). 

o 

. 1. Introduction. Random walks have been studied in abstract settings such as 

£N| ' integer lattices Z d or Riemannian manifolds ([4], [7], [10]). In applied settings there 

are many spatially explicit individual-based models (IBMs) in which the behavior of 
the system is determined by the meeting of randomly moving agents. The transmis- 
. sion of a pathogenic agent, the spread of a rumor, or the sharing of some property 

when randomly moving particles meet are examples that come to mind in biology, 
sociology, or physics ([3], [8], [6], [5], [2]). In many of these models the movement of 
agents is conceptualized as discrete transitions between square or hexagonal cells ([3]). 
However, such a stylized representation of individual movements may not always be 
entirely realistic. 

Although IBMs are powerful tools for the description of complex systems, they 
suffer from a shortage of analytical results. For example, if a susceptible and an 
infective agent move randomly in some bounded space, what is the probability of 
them meeting, and hence of the transmission of the infection? What is the average 
time until the meeting takes place? 

In the present paper we begin to answer these questions by considering a random 
walk in a bounded region of the plane or on the sphere, which we denote by R. The 
model evolves in discrete time. At each time-step an agent leaves its current position 
. at a uniformly distributed angle in the (0, 2tt) interval. On the plane the agent moves 

a fixed distance d in a straight line. On a sphere the agent moves a fixed distance d 
£\\ ■ on a geodesic. 

Here we will consider two such agents who move different distances d\ and d 2 at 
each time-step. We assume that the agents' initial positions are uniformly distributed 
in R. The paper's central results concern the probability that the paths of the two 
agents cross during the first time-step. If R is either a sufficiently large bounded 
region of the plane or a sphere, this "first-step" probability of intersection is close to 
•jjj ■ 2d 1 d 2 /(TrA[R\) where A[R] is the area of R. 

In applied settings we are often interested in the long-run average rate at which the 
paths cross. In order to extend results on the "first-step" probability of intersection 
and apply the law of large numbers we would need the following assumptions: 

• The positions of the two agents are uniformly distributed at every time step 
(which is the case on the sphere but not on the plane because of reflection 
problems at the boundary of R) , 

• The crossing-path events are independent over time (which is the case in 
neither setting because of the strong spatial dependence at consecutive time- 
steps). 
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Fig. 1.1. Example of two trajectories V\(k) and V2(k) with paths that crossed synchronously 
at mi for k = 2. (We are not interested in the asynchronous crossing at iri2 which took place at 
different times for the two agents). 

Numerous simulations have shown that despite marked departures from these 
assumptions the long-run rate at which the paths cross is also close to 2d\d,2 / '(it A[R]) . 

Section [2] contains the results both for the plane and the sphere. Section [3] is 
devoted to the numerical simulations. Extensions are discussed in Section [3J Three 
technical appendices can be found in Section [3 

2. Model. 

2.1. Geometric description in the plane and on the sphere. A bounded 
region of the plane is the most natural setting for agents moving in a 2D environment. 
However we then need to specify how agents are reflected when they hit the boundary 
of the region. There is no such problem on a sphere. In what follows R is either a 
bounded region of the plane or a sphere. 

The initial positions Vi(0) and V^O) of the two agents are assumed uniformly 
distributed in R. At each time-step the two agents move distances e?i and d-2 (which 
are fixed positive parameters) in a straight line (or along a geodesic on a sphere) . They 
depart at random angles a.\ and oli that are uniformly and independently distributed 
over (0, 2tt). The endpoints after the k — th time-step are V\(k) and V^fc) (Figure 

ED. 

The definition of a meeting in such a model is tricky because the probability of 
the two agents being in exactly the same position at any given period is 0. There are 
however different ways of approximating such a meeting. 

One could say that the agents meet if the distance between two points V\{k) and 
Va(fc) is less than some e. In such a definition results would depend on e, which is 
undesirable. For this reason we choose to define a meeting during the k — th time- 
step when paths cross between k and k + 1. This means that the segments (or the 
"geodesic arcs") Vx{k)V\{k + 1) and V2(k)V2(k+ 1) intersect. We recognize that Vi(fc) 
and V^fc) can then be close without the paths crossing, but at least the definition 
does not depend on some arbitrary e. 

The point mi in Figure [TTT1 is an example of such a synchronous crossing of paths. 
Of course the agents are not at mi at the same time. In the figure the paths cross 
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asynchronously at m-x- 

Much will depend on whether V±(k) and V 2 (k) are uniformly distributed on R for 
every k. This will be the case if R is a sphere because Vi(0) and V 2 (0) are themselves 
uniformly distributed. If on the other hand R is a bounded region of the plane, then 
the uniformity in the distributions of Vi(fc) and V 2 (k) is compromised for k > by 
the vexing problem of the behavior of the agents when they hit the boundary of R. 

For this reason we focus on the probability of paths crossing at the first time-step 

def def 

only. We simplify notations by letting V± —' (xi,yi) and V2 =' {x 2 , 2/2) be the initial 
positions of the two agents. The corresponding endpoints are denoted W\ and W 2 . 
We next proceed with calculations when R is a bounded region of the plane. 

2.2. Crossing-path probability in a bounded region of the plane. In the 

plane the endpoints W\ and W2 are 

W x d = Vx + di(cos{ai), sm(ai)) (2.1) 

W 2 d = V 2 + d 2 (cos{a 2 ), sin(a 2 j) (2.2) 

where ct\ and a 2 are polar angles that are uniformly distributed in (0, 2ir). 

We will finesse the reflection problem at the boundary of R by considering the 
possibility of intersection on the basis of Eqs. (|2.ip - (|2.2p even if W\ or W 2 is outside 
R. In such a case we examine first whether the intersection has occurred. We then 
move, in some unspecified manner, the wayward point(s) back inside R. 

We now define the feasible domain FD(Vi, a±) as the set of points V that are 
within d 2 of the segment V\W\: 

FD{Vi,ax) d = {V = (x,y) such that d(V, (^Wi)) < d 2 } (2.3) 

where d(V, (V1W1)) is the distance between a point V and the segment V\W\. 

The point V 2 must be in FD(Vi,ai) in order for the segments V\W\ and V2W2 
to intersect - although the condition is not sufficient. Figure |2~T1 depicts a segment 
V\W\ and the corresponding feasible domain FD(Vi,ai). This domain is bounded 
by dotted lines (a rectangle with sides d\ and 2d 2 with a half circle of radius d 2 at 
each end of the rectangle). In the figure a point V 2 is in the feasible domain and the 
two segments V\Wi and V 2 W 2 intersect. Whether there is an intersection depends on 
the angle a 2 . 

The region R is partitioned into an inner region Ri and a border region Rb char- 
acterized by a distance to the outside world that is either larger (for Ri) or smaller 
(for Rb) than d\ + d 2 . (The point is that FD(Vi, a±) is entirely in R if V% & Ri). We 
defined [S] as the area of a bounded subset S of M 2 . 

With V\ = (xi,yi) and V 2 = {x 2 ,y 2 ) uniformly and independently distributed on 
R, we will now calculate the probability that V1W1 and V 2 W 2 intersect. This first-step 
probability of intersection depends only on d\,d 2 and R, and is noted P p {d\, d 2 , R). 
(The subscript p indicates a probability in the plane). 

We first need the probability G p {x\ 1 yi,a\) of the intersection conditionally on 
V\ = (a;i,yi) and ol\\ G p (xi,yi,ai) is the probability (denoted by p\ below) that V 2 
falls in FD((xi,yi),ai)f]R multiplied by the probability (denoted by p 2 below) of 
the segments intersecting conditionally on V 2 falling in FD((xi, yi),oti) f] R. (In the 
F D function we have replaced V± by its coordinates (x\,y\)). 
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Fig. 2.1. Example of two intersecting segments ViWi and V2W2 with feasible domain 
FD(Vi, Inner and border regions Ri and Ri, are at distances to the outside world that are 

larger than and less than the sum di + d2- (The distance di + d^ between the two boundaries is not 
to scale). 



Given the uniformity assumption on V2, the probability p\ is then 



Pi 



A[FP((a;i,yi),ai)r|fl] 
A[R] 



(2.4) 



In order to calculate P2 we need to define the probability f p (xi,yi, a±, X2, 2/2) that 
V2W2 intersects V\W\, conditionally on V\ = {x\,yi),V2 = (2:2,2/2) and the angle a\. 
(This probability, derived in Appendix A, is obtained by calculating the magnitude of 
the angle /3 within which 0.2 must fall for the intersection to occur, and then dividing 
by 2tt) . The conditional probability p 2 of the intersection is then 



A[FD{{ Xl , yi ),a{)^R] 



The probability G p {x\,y\,a\) is now the product p\P2 which simplifies to 



(2.5) 



, , ff( X2 ,y 2 )eFD(( Xim) , ai ) nfl yi. »i. V2)dx 2 dy2 
G p (xi,yi,ai) . (2.b) 



When Vi = (xi,yi) is in Ri (i.e. the feasible domain is entirely in R) then the 
double integral on the right-hand side of Eq. (|2.6|1 is independent of V\ and is noted 
I p (di,d 2 ), i.e. 



I p (di,d 2 ) = ' // f P {xi,yi,a 1 ,x 2 ,y 2 )dx 2 dy 2 . (2.7) 

The quantity I p (di,d 2 ) is an upper bound for the double integral in Eq. (|2.6[) when 
Vi = (xx,yx) is in the border region Rb (because the integration is then over an area 
smaller than the feasible domain FD(xi, yi)). 

We now have the following result on I p (di, d 2 ), G p (xi,yi, a±), and P p (di, d 2 , R). 

Proposition 2.1. We have 

Ip{dx,d 2 )=2d x d 2 /-K (2.8) 

and when V\ — [x\, yi) is in Ri, 

G p (xi,yi,ai) = ^2[Ry ^ 2 ' 9 ^ 

When V\ = (xi,yi) is in Rb then 

G p (x liyi ,ai)<^^. (2.10) 

The first-step probability of intersection P p {d\, d 2 , R) satisfies 

def. 2d 1 d 2 A[R l ] def. 2did 2 

which leads to the mid-point approximation 

pm a j>\~ * de f- J j A[Rj]+A[R] 

P p (di,d 2 ,R) ~p = did 2 — nA ^ 2 • ( 2 - 12 ) 

The absolute value of the maximum percentage error (AVMPE) made with the 
approximation of is 

■™ = ™f™ < 2 ' 13 > 

Proof. See Appendix A. □ 

Remark. The approximation of Eq. (|2.12p is of interest and the error of Eq. (|2.13p 
is small only if A[R] is large enough in the sense that A[Ri\ is relatively close to A[R}. 
This means there is a large subset of R within which the feasible domain FD(Vi, a±) 
is entirely in R. Suppose for example that R is a circle of radius r > d± + d 2 , then 

I _ (l _ d 1 +d 2 \ 2 

AVMPE = 100 i r_) f 2 .\A) 

1+(1-^) 

which is approximately 100(<ii +d 2 )/r when d\ +d 2 is much smaller than r. Therefore 
if d\ + d 2 is one percent of the radius then the maximum error made with the estimate 
of Eq. (|2.12[) is also approximately one percent. 

We next turn our attention to the situation in which R is a sphere. 



2.3. Crossing-path probability on a sphere. When the domain R is a sphere 
of radius p we use the spherical system of coordinates where a point V is defined by 
the triplet (ft, 0, fa) of radial, azimuthal, and zenithal coordinates. 

def. 

Given initial points Vk = (p, Qk, fa), (k — 1, 2) each endpoint Wk is on the circle 
at a geodesic distance dk from Vk . The position of Wk on the circle is determined by 
an angle «fc uniformly distributed in (0, 2ir). See Appendix B for the exact derivation 
of the endpoints Wk- 

We let P s (di,d2,R) be the probability that the arcs V\W\ and V%W% intersect. 
Because the initial points Vk are uniformly distributed, the endpoints will also be 
uniformly distributed on the sphere. Therefore P s (d\ , di , R) is the probability of 
paths crossing at every time-step. 

In order to calculate P s (di ,d2,R) we proceed as before except that there is no 
border area and the double integrals are calculated in spherical coordinates. The 
feasible domain is denoted FD(0i, fa) and is defined with geodesic distances. The 
differential area element is now p 2 sin((f))d9d(f). 

We let f s (&i,4>i,ati, 62, (fa) be the probability of intersection conditionally on the 
arc V1W1 (defined by 6%, fa and ax) and on V2 (defined by 62 and fa). We also define 
the double integral 



I s (d 1 ,d 2 ,p) d = [[ f s {e 1 ,fa,a u e 2 ,fa)p 2 sin{fa)d92dfa. (2.15) 

J J(02,<p2)£FD((9 1 ,<p 1 ),a 1 ) 

The probability of the intersection conditionally on V\ = (p,8i,fa) and a± is 
independent of p, 9\ , fa and is denoted G s (d\ , e?2 , p) . We then have 

G s (d u d2, P ) = Isid ^ p) (2.16) 

where the area A[R] of the sphere appearing in the denominator is now 4irp 2 . 
The probability of paths crossing at each time-step is then 

SSJ i9i4i)eRaim2v) G s {e 1 ,fa,a 1 )^sin(fa)de 1 dfada 1 
P a (ai,d2,Jti) = 



2nA[R] 



^ I s (d 1 ,d 2 ,p) _ I s (d 1 ,d 2 ,p) 

A[R] ±irp 2 ' [ ■ ' 

When p — ► 00 the integral I s (di, (I2, p) on the sphere approaches the corresponding 
integral I p (di, ^2) = 2^1^2/71" on the plane (Eq. (|2.8|) ). 

In the next proposition we show numerically that I s (di,d2, p) is extremely close 
to 2d\d2/'K even when p is not particularly large compared to d\ and c?2- We will just 
assume that 

^<^ + ^ (2.18) 
p 2 2p p 



which insures that the feasible domain does not "wrap around" the sphere. 

6 



Proposition 2.2. With 12.18\) we have 



I 8 (di,d 2 ,p)^2did 2 /TT (2.19) 
and the probability of paths crossing at each time-step is 

(2 ,o, 

The differences between both sides of \2. 1 9\) and h2.20)) are so small that they are within 
the margins of error when calculating I s (di, d 2 , p) (and P s (d±, d 2 , R)) numerically. 
Proof. See Appendix C. □ 

The next section is devoted to simulations aimed at assessing the quality of the 
approximations derived above. 

3. Simulations. We consider a region R that is a circle of radius r = 10. At 
each time-step the two agents move distances d% — 1 and d 2 = 0.7 respectively. If an 
agent moves outside the circle we first check whether paths have crossed. We then 
move the wayward agent to a point diametrically opposed to its current position, at 
a distance inside the circle equal to the distance between the circle and the current 
position. This algorithm could be considered a 2D version of the wrapping around 
that takes place on a sphere. The goal is to try to keep the distribution of the agents as 
uniform as possible. We need this in order for the crossing-path probability to remain 
as close as possible to the "hrst-step" probability derived under the assumption that 
initial positions are uniformly distributed. 

The bounds of (|2.11[) and the approximation of (|2.12p are now used to calculate 
the low, high and mid-point approximations Pe.otPhiiP* f° r the first-step probability 
of intersection P s (di ,d 2 ,p): 



P* = d ^ Am ^l Rb] = °- 001198 (3-2) 

which translate into a maximum error AVMPE on p* of 18.42% (Eq. (|2"7H)) ). This 
error is relatively large because the sum d\ + d 2 is 1.7, which is not particularly small 
compared to the radius 10 of the circle. 

We let F(k) be the random variable equal to the average crossing path frequency 
over the first k time-steps. If the probability of paths crossing at every time-step 
were pi Q (or phi) and if paths crossing were independent events (which they are not) 
then the Central Limit Theorem would insure that for large k the corresponding 
"hypothetical frequency" -FK(fc) would be approximately normally distributed with 
meanp^o (or p hi ) and standard deviation y/ptJl - Ptoj/k (or y/phi{l —phij/k). 

Figure ETT1 depicts a simulated trajectory of the frequency F(k). We also plotted 
the hypothetical low and high intervals within which each Fh{k) would fall with 
probability 0.95. These bands shed light on the expected fluctuations of the frequency, 
in the case of independent events taking place with probability pi Q (or phi). 

This and other simulations suggest that p* gives at least an idea of the (long-run) 
crossing-path probability. 
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Fig. 3.1. Average crossing-path freguency F(k) (over the first k time-steps) for two random 
walks in a circle of radius r = 10 with d\ = 1 and d^ = 0.7 (15,000 time-steps). At each time-step 
the intervals within which the "hypothetical frequency" Ff l (k) falls with probability 0.95 are given 
for the low and high approximations pi a and p^i of P(d\, d^,p). The mid-point approximation p* 
is also plotted. 

There is less uncertainty on the sphere as we have shown that the one-step prob- 
ability of intersection P s (d\,d2,R) is extremely close to did 2 / (2n 2 p 2 ) . Simulations 
performed on the sphere (not shown) yield long-run crossing-path rates that are close 
to did,2 1 (2tt 2 p 2 ) even though the law of large numbers cannot be invoked because of 
spatial dependence over time. 

4. Discussion. The expression of Eq. (|2.12p for the crossing-path probability 
in the plane can be relatively crude. However it has the merit of simplicity and it 
improves if the area of the region R increases. 

On the sphere the crossing-path probability can be approximated very closely by 
the simple expression d\d 2 / (2ir 2 p 2 ) . We derived this expression on the basis of an 
analytical result for the plane (Appendix A) , expecting it to be a good approximation 
only for a sphere of infinitely large radius. It is of some interest to note that because 
of the complexity of the multiple integrals in spherical coordinates we saw now way of 
obtaining this approximation from the calculations performed directly on the sphere 
(Appendix C). 

Because of marked departures from required assumptions the law of large numbers 
could be applied neither in the plane nor on the sphere. Despite that, the long-run 
average crossing-path probabilities appear to be close to the first-step probabilities. 
This suggests that a weaker version of the law of large numbers may be applicable. 
For example results on "weakly dependent" random variables (i.e. variables that are 
"m" (or " (/?" )-dependent ([!])) may provide more insights into the long-run behavior 
of the system. 

We note some obvious and some less obvious extensions that can be of use to 
population biologists (and perhaps others): 
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Fig. 5.1. Intersecting segments ViWj and V2W2 in coordinate system in which V1W1 lies on 
the x-axis and the origin is at the middle ofV\W\. 

• If the crossing of paths takes place between an infected and a susceptible 
agent, then the transmission of the infection may occur with only a probability 
t. In such a case the crossing-path probabilities found here need simply be 
multiplied by r in order to obtain the probability of transmission at each 
time-step. 

• An important extension would have I such infectives and S such susceptibles. 
Epidemiologists would be keenly interested in analytical results on the rate 
at which the infection would then spread. 

• The assumption that an agent moves at an angle that is uniformly distributed 
in (0,27r) may not be realistic. For example animals may move only within 
a limited angle in the continuation of the previous direction. Preliminary 
investigations suggest that the results obtained here may still be applicable. 

The results given in this paper are merely starting points for more in-depth theo- 
retical investigations. They also provide practitioners with some answers concerning 
the dynamics of a process that depends on randomly moving agents meeting in a 
spatially explicit environment. 

5. Appendices. 

5.1. Appendix A: Proof of Proposition I2TT1 The double integral I p {d\,d2) 
in Eq. (|2.6p is calculated in the orthonormal coordinate system (a;, y) for which the 
segment V{W\ lies on the x-axis and the origin is at the middle of the segment (Figure 

EH) . 

In the new coordinate system, the probability fp{x\,yi, ai,%2, 7/2) that V{Wi and 
V2W2 intersect depends only on the components (£2,2/2) of V2. If this probability is 
denoted F p (x2,y2), then I p (di,d 2 ) is 



py=d 2 rx=d 1 l'l+yj d\-y 2 

I p (dx,d 2 ) = / / F p (x,y)dxdy. (5.1) 

J y— — di J x — — di/2 — yd§— : y^ 

The probability F p {x2, 2/2) is obtained by calculating the magnitude of the angle 
P within which V2W2 (defined by the angle 0.2) must fall and then dividing by 2-7T (see 
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Figure l5~Tj) . In the Figure the point V 2 is in neither circle and the angle f3 is the vertex 
angle in the isosceles triangle with apex V 2 and two sides of length d^. Therefore, 
when V 2 is in neither circle the probability F p {x 2 ,y 2 ) of an intersection is given by 
the function 

o f\yz\ 

2 arccos I — — 

F P A^y 2 ) d ^ A^-. (5.2) 

Similar geometric considerations show that if the two circles overlap (i.e. d\/2 < d 2 ) 
then for a point (x 2 ,y 2 ) in the intersection of the two circles, the probability of 
intersection F(x 2 , J/2) is given by the function 

t , x 2 +d 1 /2\ fx 2 -di/2 
arctan ; — ; — arctan ' 



F p , 2 {x 2 ,y 2 ) = —z- " -• ( 5 - 3 ) 

Finally, if (x2,y2) is in only one of the circles, then one of the vertices of the 
triangle with apex V 2 will be the center V\ or W\ of that circle. The probability of 
intersection F{x 2 ^y2) is then 

arccos f _ arctan (fl^ 



Long but elementary calculations show that the innermost integral in (|5.ip . con- 
sidered a function H (y) of y, is now equal to 

d i x arccos ( ¥ 

H(y) *M- / F p (x,y)dx = (5.5) 

Jx = -d 1 /2-y / dl-y 2 T 

We therefore have 

rV=d 2 nj j 

I p (d 1) d a )= / H(y)dy==^ (5.6) 

Jy=-d 2 TT 



which is Eq. (|2.8p and yields Eq. 

The inequality in (|2.10[) results from the fact that I p (di,d2) is an upper bound 
for the double integral in Eq. (|2.6|) when V\ — (xi, y\) is in the border region R\,. 

With V\ and V 2 uniformly distributed on R, the probability of intersection 
Pp(dx, d 2 , R) is now obtained by integrating G p {x\, y±, a±) over (x±, y±) in R — Ri [J Rb 
and over a± in (0, 2tt), and then dividing by 2nA[R\: 

, m ///(xx, yi )efl,aie(o,!hr) Gr i'( a! i'I'i> a i)dafidl'id«i 
P P (d u d 2) R) 



///(^,m)gi?„a 1 g(o,2.r) G P (a:i,yi,ai)da;idyidai 



27rA[i?] 
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2ttA[R] 



_ 2d 1 d 2 A[R l ] nS {xi , Vl ) & R b , a ^{ ^) G p^y^ a ^ dx ^ da ^ , , 
nA[R} 2 2ttA[R] ' 1 ' j 

The first term is easily calculated while the second one is a complicated triple 
integral of a double integral that can only be calculated numerically. The inequality 
in (|2.10p implies however that 

2d 1 d 2 A[R l ] 2d 1 d 2 (A[R i ] + A[R b \) 2d x d 2 
irA[R} 2 ~ P ^ d ^ R > ^ ^A[Rj2 = ^WY ( ' 

which is (|2.11|) . The approximation of (|2.12p is simply the mid-point of the interval 
in (|5.8[) . The error term of Eq. (|2.14p is based on the low and high bounds of ([57 



5.2. Appendix B: Derivation of Wi,W 2 . In order to find each endpoint Wk 
{k = 1, 2) we will first determine a particular point Qk at a geodesic distance dk from 
Vfc. Each Wk is then obtained by rotating Qk by a uniformly distributed angle cxk 
about the OVk axis. 

Each Qk is defined as having the same azimuthal coordinate 9k as V& = (p, 9k, 4>k) 
and a zenithal coordinate <p that puts Qk at a geodesic distance dk from Vfc. To 
calculate Qk we define the function 



z(<M) d = I t + d J, P [ [t <n - d /* (5.9) 
[ (p—d/p otherwise. 

Then 

Qk d = (p,O k ,z(<l>k,d k )). (5.10) 

To obtain Wk from Qk we define the normalized vectors Uk d =' OVk/\\OVk\\ with 
Cartesian coordinates Uk,x,Uk,y,Uk. z - We next define the antisymmetric matrix 



def ( _U M u k,y \ 
A k = U k ,z -U k ,x ■ (5-11) 

\ -Uk, v Uk, x / 

We let I3 denote the 3x3 identity matrix. We will also need the Cart(p, 9, (f>) and 
Sph{x, y, z) functions that transform spherical coordinates into Cartesian ones, and 
vice-versa, i.e. 

Cart(p, 9, 4>)) = p (cos(0) sin(0), sin(0) sin(0), cos(0)) (5.12) 

and 



Sph(x,y,z) = f V x 2 + y 2 + z 2 , 9{x, y, z), arccos ( ^ - ^ 2 - ^ J J (5.13) 
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where 

, . | arcsin , y if < x 

qu y z \ d g- ) V* 2 +v 2 (5 14) 

I 7r — arcsin , y otherwise. 

Each Wk, expressed in spherical coordinates, is now obtained by using Rodrigues' 
rotation formula, (pQ) i.e. 



W k Sph ([I 3 + A k sin(a fe ) + A 2 k (l - cos{a k ))]Cart(Q k )) (5.15) 

where the a k 's are uniformly distributed in (0,27r). 

5.3. Appendix C: Proof of Proposition 12.21 We will derive an expression 
for I s (di , d,2 , p) in the spherical coordinate system in which the arc V\Wi lies on the 
equator and its middle is at the point (p, 0, 7r/2), i.e. 

Vx = (p, -d 1 /(2p), tt/2), Wi = (p, d 1 /(2p), tt/2). (5.16) 

The feasible domain now consists of two areas on the surface of the sphere. First the 
"spherical rectangle" centered on (p, 0, tt/2) with sides of geodesic lengths di and 2d2', 
and second at both ends of the rectangle the half-circles centered at V\ and W\ and 
of (geodesic) radius di (Figure EH]). The feasible domain now depends only on d\,di 
and is denoted FD(d\,d2)- 



0.25 -, 




Fig. 5.2. Feasible domain on the sphere and angle {3 within which V2W2 must fall for the arcs 
V1W1 and V2W2 to intersect. (When V2 is in the right circle (resp. the left circle), h r (resp. hi) 
will be at W± (resp. V\). 

We need several functions in order to calculate G s (9i, 4>i,ati): 
• The geodesic distance function between two points Pi = (p, 9%,<f)i) and P2 = 
(p, 02, 4>2) on the sphere: 

gd(P 1 ,P 2 ) d ^ 
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parccos[sin(</>i) sin((/)2)(cos(0i) cos(0 2 ) + sin(#i) sin(6>2)) + cos(0i) cos^)]- 

(5.17) 

• If X and Y are two vectors (in Cartesian coordinates) on the sphere of radius 
p, then the normed tangent vector at X to the geodesic line between X and 
Y: 



t(X,Y) 



def. 



Y 



X 1 Y 

\\xv 



X 



(5.18) 



• The Cart(p, 9, cj>) function that transforms spherical coordinates into Carte- 
sian ones (Eq. (|5.12p ). 
Given V2 = (p, O2, <fa) in the feasible domain FD(dx, ^2), we let h r and hg be the 
two points on the equator (on the right and on the left of V2) that determine the 
magnitude of the angle (3 within which V2W2 must fall for the intersection to occur 
(Figure 

Bearing in mind Vy and W\ of Eq. (|5.16p . the Cartesian coordinates of h r and 
hi, considered functions of V2 — $2, $2), are 



h r {92,4>2 



def. 



Cart(Wi) if gd(V 2 , Wx) < d 2 

cos^/p) \ . / 
sm{<f> 2 ) J 1 \ 



p I cos I W 2 + arccos ■ 



2 + arccos c °^^ ) .()) otherwise 



(5.19) 



hi{d2,<fe) 



def. 



Cart{Vx) \fgd{V 2 ,Vx) < d 2 
arccos — . \r T ; 

sm(0 2 ) 



(cos (6 



^ , sin (b 



cos(rf 2 /p) 
sin(0 2 ) 



] otherwise. 



(5.20) 



In the new coordinate system, the probability of intersection / s (#i, <f>i, cn\, 92, </> 2 ) 
conditionally on the arc V{W\ and on V2 depends only on the components {p, 92,4>2) 
of V2 and is denoted F s {92,4>2)- Given the uniformity assumption for the angle a.2 
at which the second agent leaves V2 to go to W2, the probability is equal to the 
magnitude of the angle j3 within which the arc V2W2 must fall for the intersection to 
occur, divided by 2tt. The angle (3 is the angle between the tangent vectors at V2 in 
the directions of ^(#2, </> 2 ) and h r (92,4>2)- The probability of intersection -F s (# 2 ,0 2 ) 
is therefore 



arccos (r[Cart(p, 9 2 , 2 ), h r (9 2 , (j)2)] T r[Cart{p, 6 2 ,<j>2), h e (9 2 ,(/)2)} 

"hi 



(5.21) 

The double integral I s {d\ ,d2,p) is now calculated in the new coordinate system 
by integrating F s {92, $2) over the feasible domain FD{d\, g^). Because of symmetries 
the integral I s {d\, c?2, p) is the sum of four times the integral over the upper right 
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Table 5.1 

Percentage error 100 X ( i^d^d^p) ~ when approximating I s (di,d2, p) as 2did2/ir, with 
di = 3 and an illustrative range of values for d2 and p. 



d>2 4;p-» 


2 


3 


4 


5 


1 


4.866 x icr 4 


-1.68 x 1(T 4 


-2.527 x 1CT 4 


-8.471 x 1(T 4 


2 


1.288 x 1(T 3 


5.595 x 1CT 4 


9.269 x 1CT 4 


6.658 x 1(T 4 


3 


-1.506 x 1CT 5 


7.529 x 1CT 4 


9.412 x 1CT 5 


8.954 x 10~ 4 



quarter of the rectangular area of FD{d\, d%) and of four times the integral over the 
upper half of the right circle. We thus have 

rd 1 /(2p) mt/2 

I s (d 1 ,d 2 ,p) =Ap 2 / F s (62,<j>2)sin(<fa)d6 2 d(t>2+ 

Je 2 =o J<t>2=n/2-d 2 /p 

rd 1 /(2p)+d 2 /p pn/2 

4p 2 / / F s (9 2 ,<\> 2 )sin{ch 2 )d9 2 d<f, 2 (5.22) 

Je 2 =d 1 /{2 P ) J <t> 2 =<t>{0 2 ) 

where 




(5.23) 



is the lower value of 4> 2 when integrating in the right half circle. 

Over a wide range of values for d\ , d 2 and p we found that the relative error made 
when approximating I s {d\, d 2 ,p) of Eq. 1|5.22|) as 2did 2 /ir was of the order of 10 -3 to 
10 -5 percent. See Table I5TT1 for an example of this relative error for a range of values 
of d 2 and p. In fact the error is so small that it is within the margin of error when 
calculating I s (di , d 2 , p) numerically. The fact that I s (d\,d 2 , p) ~ 2did 2 /ir combined 
with Eq. ([2~T71) yields the result of l[2"3U]) , 
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